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LONG-TERM  GOALS 

Central  to  the  long  term  goals  of  this  joint  project  is  to  understand  the  physics  of  the  propagation  of 
uncertainty  through  the  interfaces  between  oceanography,  acoustics,  array  processing  and  performance 
prediction.  We  will  develop  an  efficient  overall  simulation  platform  that  combines  all  of  the 
components  of  the  baseline  (mean)  and  uncertainty  problem:  Oceanography  through  4D  acoustic  field 
prediction.  The  development  of  a  methodology  to  distill  the  complexity  and  uncertainty  of  the  ocean 
acoustic  environment  and  the  system  level  sensitivities  to  relevant  situational  awareness  for  the 
operator  is  an  important  goal  of  this  research. 

OBJECTIVES 

The  objective  of  this  research  program  is  to  both  develop  a  methodology  to  predict  uncertainty  in  the 
whole  performance  prediction  process  and  to  understand  the  uncertainty  physics  of  the  individual 
components  of  the  process.  The  latter  provides  the  potential  to  develop  methods  to  reduce  uncertainty. 
We  intend  to  follow  a  two  pronged  approach:  (1)  total  model  development  for  Monte  Carlo  simulation 
and  (2)  studying  the  physics  of  the  interfaces  between  oceanography,  acoustics,  array  processing  and 
performance  prediction. 

APPROACH 

This  project  is  a  joint  effort  between  the  parties  listed  above. 

We  intend  to  follow  a  two  pronged  approach:  (1)  total  model  development  for  Monte  Carlo  simulation 
and  (2)  studying  the  physics  of  the  interfaces  between  oceanography,  acoustics,  array  processing  and 
performance  prediction. 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number 


1.  REPORT  DATE 

30  SEP  2001 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2001  to  00-00-2001 


4.  TITLE  AND  SUBTITLE 

Modeling  and  Analyzing  the  Propagation  from  Environmental  through 
Sonar  Performance  Prediction 

6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Scripps  Institution  of  Oceanography „La  Jolla„CA,  92093  report  number 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Central  to  the  long  term  goals  of  this  joint  project  is  to  understand  the  physics  of  the  propagation  of 
uncertainty  through  the  interfaces  between  oceanography,  acoustics,  array  processing  and  performance 
prediction.  We  will  develop  an  efficient  overall  simulation  platform  that  combines  all  of  the  components  of 
the  baseline  (mean)  and  uncertainty  problem:  Oceanography  through  4D  acoustic  field  prediction.  The 
development  of  a  methodology  to  distill  the  complexity  and  uncertainty  of  the  ocean  acoustic  environment 
and  the  system  level  sensitivities  to  relevant  situational  awareness  for  the  operator  is  an  important  goal  of 
this  research. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a  REPORT 

unclassified 

b  ABSTRACT 

unclassified 

c  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

5 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


(1)  [Cornuelle,  Kuperman,  Hodgkiss,  Porter/Hursky,  Cox/Heaney]  Our  primary  approach  to  the  total 
uncertainty  simulation  model  development  will  depend  on  the  usual  Monte  Carlo  runs  to  convert  an 
ensemble  of  oceanographic  states  to  an  ensemble  of  system  output  showing  the  range  of  possible 
values.  This  method  has  the  advantage  of  retaining  validity  in  the  presence  of  strong  nonlinearity,  as 
well  as  being  computationally  simple. 

The  ocean  ensemble  will  be  coupled  to  a  wide  area  acoustic  propagation  model  and  an  accurate  array 
processing  model  for  the  purpose  of  perfonning  area  wide  perfonnance  prediction  with  an  uncertainty 
measure.  Particular  care  will  be  given  to  the  mode  of  representing  this  information  to  the  operator. 

We  also  intend  to  explore  linearized  methods  for  determining  system  response  ranges,  based  on  an 
analogy  to  the  tangent  linear  model  and  adjoint  model  techniques  in  use  in  Physical  Oceanography. 
Applying  these  concepts  (and  possibly  something  like  the  automated  Giering  compiler)  to  the  acoustic 
codes  may  provide  a  computationally  efficient  way  to  compute  sensitivities  and  transform  ocean  state 
uncertainty  covariances  to  system  perfonnance  uncertainty  covariances.  It  may  be  that  only  an 
idealized  reduction  of  the  problem  (as  set-up  in  1)  is  attempted,  depending  on  the  difficulties  that  are 
uncovered.  Further,  it  is  hypothesized  that  this  approach  may  be  more  applicable  to  the  total  problem 
rather  than  stopping  at  the  acoustics  output  because,  to  a  certain  extent,  performance  prediction  is 
somewhat  of  a  “smoothing”  process,  possibly  enabling  the  required  linearization  mentioned  above.  We 
consider  this  part  to  be  a  high-risk  high-payoff  idea,  in  that  to  our  knowledge,  the  oceanographic  data 
assimilation  adjoint  approach  has  never  been  applied  to  acoustic  models,  much  less  to  the  whole 
performance  prediction  process. 

Uncertainty  occurs  in  the  geometry  (source/receiver  locations,  bottom  depth),  the  bottom  properties, 
the  surface  properties  (sea  state  and  resulting  bubble  clouds  and  surface  roughness),  and  the  ocean 
volume  (internal  waves/tides,  meso-scale  features  such  as  eddies,  and  fronts).  The  uncertainty  of  many 
of  these  is  easy  to  characterize.  For  instance  errors  in  receiver  locations  are  a  function  of  the 
performance  of  depth-heading  sensors  on  the  arrays;  their  perfonnance  may  vary  over  time  on  a  given 
array  but  by  and  large  there  are  no  interesting  scientific  issues  about  their  accuracy.  The  most 
interesting  and  perhaps  important  of  the  above  listed  uncertainties,  is  that  due  to  the  oceanography  and 
how  it  subsequently  regulates  acoustic  bottom  interaction  with  uncertain  geophysical  parameters. 

The  ROMS  primitive  equation  model  will  be  used  to  simulate  the  small-scale  features  in  an  area  of 
interest.  Separately,  a  simpler  internal  wave  model  will  be  used  to  characterize  the  space-time 
structure  of  the  internal  wave  field  that  will  be  excited  using  a  unifonn  energy  density  of  internal 
waves.  In  an  operational  scenario,  we  envision  that  one  or  the  other  of  these  models  will  be  used  to 
predict  the  statistics  of  the  oceanographic  variation.  The  primitive  equation  model  is  not  viewed  as  a 
tool  that  would  predict  a  deterministic  environment  but  rather  as  something  that  might  provide  a  more 
accurate  prediction  of  the  typical  variation.  The  1-D  internal  wave  model  is  obviously  simpler  and 
computationally  more  practical,  but  may  be  less  accurate.  One  objective  will  be  to  compare  the 
relative  merits  of  the  two  approaches. 

The  existing  acoustic  models  produce  a  realization  of  the  pressure  field  for  a  single  deterministic 
environment  (which  is  obviously  distinct  from  a  mean  field).  To  capture  the  uncertainty  the  acoustic 
models  will  need  to  be  enhanced  to  rapidly  produce  an  ensemble  of  pressure  fields  or  statistics  of  the 
ensemble.  The  algorithmic  approach  we  will  develop  in  this  program  is  different  in  each  of  the  4 
standard  model  types.  However,  the  common  starting  point  for  all  of  these  is  the  ’environmental 
endpoints',  i.e.  the  limits  (or,  more  precisely,  variances)  characterizing  the  uncertain  environment. 
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Despite  the  variety  of  sources  of  uncertainty,  they  can  all  be  treated  using  the  same  framework  as  will 
be  discussed  below.  To  fix  ideas,  imagine  a  mean  sound-speed  profile  and  a  lowest-order  EOF 
characterizing  the  variation  due  to  the  first  baroclinic  mode.  The  ’environmental  endpoints'  are  the 
mean  with  that  EOF  added  and  subtracted  based  on  the  excursion  seen  in  the  oceanographic  data.  If 
the  internal  wave  model  is  found  to  be  adequate,  it  will  be  integrated  with  the  acoustic  model  so  that 
an  input  SSP  is  used  first  to  drive  the  internal  wave  model,  and  then  passed  directly  to  the  acoustics 
model  along  with  the  environmental  endpoints  calculated  by  the  internal  wave  model. 

We  also  expect  to  use  the  tools  that  we  have  developed  to  study  the  viability  of:  (a)  assimilation  of 
acoustic  data  from  ships  of  opportunity;  this  might  involve  data  fusion  with  remote  sensing  or  the  use 
of  battlegroup  location  via  netcentric  operations  and  (b)  optimizing  ASW  area  coverage  in  a  way  that 
reduces  uncertainty  by  accumulating  and  assimilating  data  during  operations. 

(2)  [Krolik]  The  results  from  the  total  simulation  approach  presented  above  will  be  difficult  to  interpret 
in  terms  of  the  individual  physics  components  of  our  program.  Hence,  there  will  be  a  need  to  study  the 
interfaces  between  the  components.  We  will  develop  computationally  efficient  mappings  of  random 
oceanographic  state  variables  into  representations  of  the  acoustic  pressure  field  received  at  a  sensor 
array.  The  goal  is  to  accurately  capture  the  spatial  wavefront  uncertainty  of  submerged  and  surface 
acoustic  sources  in  order  to  improve  in  situ  passive  sonar  performance  prediction.  The  emphasis, 
therefore,  will  be  to  develop  computationally  fast  methods  for  mapping  high-dimensional  random 
oceanographic  state  variables  to  functions  of  the  acoustic  wavefront  used  to  form  target  detection 
statistics.  Both  the  classical  complex  multivariate  Gaussian  wavefront  model  and  the  proposed  non- 
Gaussian  pseudo-multipath  expansion  (PME)  model  will  be  evaluated  for  different  passive  sonar 
scenarios  in  uncertain  oceanographic  environments.  To  benchmark  these  fast  solutions,  sonar 
performance  predictions  for  limited  complexity  cases  will  be  compared  to  those  obtained  by  Loren 
Nolte  (see  below)  using  Monte-Carlo  estimation  techniques.  Oceanographic  uncertainty  will  be 
characterized  by  the  mean  and  covariance  of  mesoscale  state  variables  generated  by  Comuelle’s 
extensions  to  the  Regional  Ocean  Modeling  System  (ROMS)  and  statistical  characterizations  of 
internal  wave  variability  developed  by  Heaney  Computational  models  developed  by  Porter  and 
Kuperman  will  be  used  to  generate  realizations  of  the  acoustic  field  in  highly  inhomogeneous  range 
and  azimuthally  dependent  littoral  environments. 

[Nolte]  The  acoustics/signal  processing  interface  involves  understanding  the  optimal  detection, 
classification,  and  localization  framework  that  incorporates  uncertainties  in  ocean  environmental 
parameters  directly  into  the  computation  of  improved  sonar  perfonnance  prediction.  Beyond  merely 
predicting  performance  degradation  due  to  environmental  uncertainties,  we  also  seek  to  identify  those 
physical  parameters  whose  uncertainties  affect  detection  performance  prediction  the  most.  By 
formulating  hypothesis  testing  problems  which  properly  incorporate  the  effects  of  uncertain 
parameters,  we  will  provide  performance  measures,  such  as  Receiver  Operating  Curves  (ROC's)  and 
localization  ROC's,  which  more  accurately  include  the  uncertainties  of  the  realistic  detection  problem 
faced  by  an  operator  in  a  complex  littoral  environment,  than  does  the  classic  sonar  equation.  The  input 
to  this  fonnalism,  is  mesoscale  and  internal  wave  simulations  from  above  as  well  as  the  acoustic 
propagation  simulations  as  passed  through  the  Krolik’ s  wavefront  distortion  representation.  This 
“data”  will  be  incorporated  into  the  fonnulation  of  uncertainties  for  our  hypothesis  testing  approach  to 
sonar  performance  prediction.  This  will  utilize  and  complement  the  statistical  characterization  of  the 
uncertainty  of  physical  parameters,  per  se,  as  determined  by  the  ocean  modelers.  In  addition  to  IMAT 
like  outputs,  various  presentations  to  the  operator  will  be  considered,  such  as  probability  of  detection 
contours  versus  range/dep th/bearing,  for  fixed  probability  of  false  alann,  as  a  function  of  ocean 
environmental  uncertainties. 
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Throughout  the  program,  we  will  have  as  underlying  theme  and  goals,  the  construction  of  a 
methodology  for  capturing  uncertainty  in  perfonnance  prediction  from  the  overall  system  perspective. 
The  focus  of  this  work  will  be  distilling  the  complexity  and  uncertainty  of  the  ocean  acoustic 
environment  and  the  system  level  sensitivities  to  relevant  situational  awareness  for  the  operator.  [Cox 
will  lead/coordinate  this  aspect  of  the  program]. 

WORK  COMPLETED 

This  project  is  just  beginning. 

RESULTS 

This  project  is  just  beginning. 

IMPACT  /  APPLICATIONS 

This  expected  impact  of  this  project  is  to  provide  a  methodology  to  provide  a  reliability  measure  to  the 
operator  of  at-sea  performance  prediction  models. 

TRANSITIONS 

No  transitions  took  place  in  FY01. 

RELATED  PROJECTS 

This  is  one  of  the  programs  in  the  ONR  UNCERTAINTY  DRI. 

PUBLICATIONS 

None 
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